GPS test of the local position invariance of Planck's constant 



(N 

o 



7'- 



> 

O 



en 

o 

(N 



X 



J. KentoslE and M. Mohagejll 

Department of Physics and Astronomy, California State University, 

Northndge, Northndge, California 91330-8268, USA 

(Dated: February 29, 2012) 

Publicly available clock correction data from the Global Positioning System was analyzed and 
used in combination with the results of terrestrial clock comparison experiments to confirm the 
local position invariance (LPI) of Planck's constant within the context of general relativity. The 
results indicate that h is invariant within a limit of |/?h| < 0.007, where Ph is a dimensionless 
parameter that represents the extent of LPI violation. 
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Many experiments and observations have tested possi- 
ble variations of fundamental constants, particularly the 
fine-structure constant a jl|. Though most results have 
been null, studies of quasar absorption spectra suggest a 
small change in a over the past 8 billion years across the 
observable universe [3l. If a were to vary as suggested 
by those results, that would imply a change in at least 
one of the parameters that comprise a = e'^ /Aireofic 3]. 
In this letter, data from the Global Positioning System 
(GPS) is analyzed in combination with prior terrestrial 
clock comparison experiments 4] to establish bounds on 
the invariance of Planck's constant h. 

This analysis provides a test of 'local position invari- 
ance' (LPI), which states that in local, freely falling 
frames the outcome of any nongravitational experiment 
is independent of where and when in the universe it is 
performed [5| . This principle is a variation of Einstein's 
equivalence principle relating gravitational fields to the 
accelerations of bodies in free space. Nevertheless, as a 
macroscopic theory of gravity, general relativity tells us 
little about h. Thus, experimental evidence of its invari- 
ance is necessary. 

The Standard Model Extension (SME) has been devel- 
oped partly to study possible variations in fundamental 
constants [l|, la, LZf . A less broad approach is used herein, 
based on the assumption that general relativity correctly 
describes macroscopic phenomena. The observable ef- 
fects of a variable h are deduced from the invariance of 
the proper energy of atomic transitions. This approach 
provides a consistency check of general relativity, and our 
results may be relevant to the review of recent observa- 
tions of super-luminal neutrinos [8| . 

The analysis proceeds in three steps: First, we use GPS 
data to determine LPI limits for atomic clocks. Then we 
incorporate the results of clock experiments to determine 
LPI limits for time dilation. Finally we use both results 
to estimate the invariance of Planck's constant. 

Using the standard convention (e.g., [9|) for expressing 
violations of LPI to first order in a gravitational field, 
the atomic transition frequency used in an atomic clock 
at X and at some arbitrary reference point O - such as a 



terrestrial laboratory - can be related as follows: 

/W/o = l + (l + /3/)AC//c2, (1) 

where fo is the frequency measured by an observer at 
point O when the clock is at O, and f^o is the frequency 
when the clock is at x as measured from the reference 
frame at O. The potential difference is AU = Ux — Uo, 
where Ui — $i — ff /2, <I>i is the gravitational potential 
energy per unit mass and Vi is the clock's velocity. The 
dimensionless parameter /3f represents the extent of LPI 
violation for atomic clock frequencies. 

The GPS uses atomic clocks on the ground and in or- 
bit. The GPS is operated on GPS time, a continuous 
timescale that can be related to Coordinated Universal 
Time (UTG). To correct for clock drift, a clock correction 
for each satellite is broadcast as part of the navigation 
signal. GPS satellite orbits have eccentricities generally 
less than 0.02 'lO[. At apogee, a satellite has a slower 
speed and higher gravitational potential, which cause its 
clock to run faster than at perigee and faster than a clock 
on the ground. Corrections for this effect are made in 
GPS, based on general relativity 10l.lll|. 

We analyzed GPS satellite data made available on the 
Internet by the International GNSS Service (IGS) [13]. 
The IGS 'Final' product in SP3-C format was used, in 
which GPS satellite positions and clock corrections are 
published for every 15 minutes of GPS time. Positions 
are precise to 1 mm and clock corrections are precise to 
1 ps [l3(. A single SP3 file includes position and clock 
correction data for 32 GPS satellites for one GPS day. 
The information in the files is based on data collected by 
IGS from the 32 satellites, 12 GPS control stations, and 
approximately 350 ground stations 
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Of the 32 GPS satellites, seven were chosen that had 
more eccentric orbits and the most stable clocks: 2, 11, 
18, 21, 26, 28, and 32. Those satellites had the lowest 
ratios of the standard deviation of their changes in clock 
corrections to the theoretical range in clock rates for that 
satellite. Satellite eccentricities ranged from 0.01 to 0.02. 

The clock corrections published by IGS do not include 
the general relativistic corrections for the eccentricities of 
the GPS satellites. It is left to the user to calculate and 



add those effects [l5| . Since we did not have to subtract 
it out, the missing relativistic component to the clock 
corrections simphfied the analysis. 

To find Pf the changes in the IGS-published clock cor- 
rections of the GPS satellites over each 15 minute inter- 
val were analyzed. One year of such data was reduced, 
spanning the time period from April, 2010 through May, 
2011, including extra days used for averaging. For each 
satellite, 95 changes in its clock correction during each 
day were plotted as a function of its distance from the 
earth's center. (The 96th spanned two SP3 files and was 
omitted for simplicity.) A representative plot is shown in 

Fig.m 
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FIG. 1: Changes in IGS-published clock corrections every 15 
minutes for Satellite No. 18 on 6/10/10, plotted against dis- 
tance from the earth's center. The uncertainty in the ordinate 
is a, as published in the SP3 files. The least-squares linear fit 
has a small slope m. 

Since the relativistic eccentricity effect is not included 
in the published clock corrections, changes to those cor- 
rections should not vary with distance from the earth's 
center. Standard theory predicts the slope m of the linear 
fit in Fig. [T]to be zero. The small slope evident in Fig. 
[I] extracted by a linear least-squares regression fit to the 
data, might be explained by a non-zero /3/. A non-zero 
slope could also be caused by any number of effects on 
satellite clock corrections, including random errors and 
atmospheric effects. 

To estimate conservatively large bounds for /3f, all 
nonzero slopes were attributed to /3f. The slope m in 
Fig. [T]is related to a possible non-zero value of f3f by 



/3/ = m{ra ~ rp)/{2ATrrMx), 



(2) 



where Tq is the apogee radius, r^ is the perigee radius 
and ATjnax is the maximum theoretical clock correction 
of a satellite at apogee, given by the following, adapted 
from 0: 

ATmax = Tep ($a " vl/2 - $GPS + ^'gPs/S) / C^ ■ (3) 

In this equation, $a and Va are the gravitational potential 



per unit mass and velocity of a satellite at apogee, respec- 
tively; and $GPS and wgps are the same for a satellite 
in an ideal, circular GPS orbit corresponding to a radius 
of 26,561.75 km, representing GPS time. (GPS satellites' 
clocks are adjusted to run more slowly pre-launch so that 
when they reach orbit, they run at GPS time.) T^p is the 
'epoch time,' which is the 15 minute interval between 
clock corrections. 

The coordinates in the SP3 files are published in an 
earth-fixed reference frame. A satellite's velocity rela- 
tive to the earth's non-rotating frame was interpolated 
from its instantaneous radial distance, using Newtonian 
gravity. For an ideal, circular GPS orbit, this method 
reproduces the standard velocity ugps of 3,873.83 m/s, 
used as a baseline in Eq. (3). 

The mean offset of the changes in clock corrections 
from zero in Fig. [1] indicates how much faster or slower 
a satellite's clock is advancing than GPS time. That has 
no bearing on the current analysis. Of interest is how 
the changes in clock corrections might vary with radial 
distance from the earth's center. 

The best fit slope and the corresponding daily values 
for (3f were calculated for the seven satellites, for each day 
in the test period, with 2,749 plots generated. (The clock 
data was partially missing or corrupted on some days.) 
Results for three representative satellites are shown in 
Fig. Ha-c). 
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FIG. 2: Daily and filtered values of /?/ for one year for (a) Sat 
2, (b) Sat 18, (c) Sat 32. Sat 2 had the least annual oscillation; 
Sat 32 had the greatest range in j3f . (d) Combined probability 
distribution of filtered /?/ for 7 satellites for one year. 

The plots of /?/ exhibit significant daily variations, and 
an annual oscillation that varied among satellites. The 
daily variations were consistent in magnitude among the 
satellites. The annual oscillations are likely caused by 
atmospheric effects. The maxima of these oscillations 
roughly coincide with a semi-major axis aligned with 



Eurasia and the Pacific Ocean, an orientation that shifts 
360 degrees over a year. 

Time-series and Fourier analysis of the data was used 
to develop appropriate digital filters to remove yearly, bi- 
monthly and daily variations driven by nongravitational 
effects. To remove the annual oscillation, a least-squares 
fit was computed using y ~ A+B cos(27ri/lyr-f-^), where 
A, B and (j) are fit parameters. The B cos{2iTt/lyr + 4>) 
term was subtracted from the daily values. Based on the 
frequency distribution in the Fourier analysis, we then 
used a 51-day moving average to filter out short-term 
fluctuations in /3/. The filtered data for three satellites 
(2, 11, 32) exhibited a sawtooth-shaped oscillation with 
a 120-day period, of unknown origin. That oscillation 
was filtered out for those three satellites in the same way 
as the annual oscillations. The limits on Pf were most 
sensitive to filtering out the annual oscillations, but less 
sensitive to the duration of the rolling average or to re- 
moving the 120-day oscillations. The residual daily val- 
ues of Pf for three satellites are shown with solid lines in 
Fig. Ha-c). 

The residual values of /3/ from seven satellites exhibit 
the bell-shaped distribution shown in Fig. WiA)- ^ small 
residual mean of /3/ = 0.0003 was present. We cannot 
rule out GPS systemic errors as its source. Based on 
a 95% confidence level, the following limits were found: 
-0.0027 < /3/ < 0.0033. The resuhs for the seven satel- 
lites were weighted equally. Finding these limits on /3/ 
is the first step in setting limits on the invariance of 
Planck's constant. 

The well-behaved annual oscillations in Fig. [2] demon- 
strate some statistical significance for the linear fits for 
771. However, values of B? for the fits ranged from 10~^ 
to 0.6 with a mean of 0.04, suggesting a low to moderate 
statistical significance for m [17 1. 

To study this effect, a statistical model was developed 
in which the standard deviations of the GPS clock correc- 
tions were used to create a randomly-generated data set 
to mimic the GPS data. Ninety-five randomly generated 
points were used to generate fits for m, repeated 2,555 
times and then filtered to model our analysis. The results 
indicate that about 2/3 of the limits on Pf arise from ran- 
dom variations in the GPS clock corrections. However, 
some GPS satellites exhibit more consistent clock correc- 
tions and the bounds could be reduced by analyzing only 
those satellites. 

If h were to vary, general relativistic time could differ 
from atomic time. Possible LPI violations of gravita- 
tional time dilation can be expressed to first order as 



dt^/dto = l + {l+Pt)^U/^ 



(4) 



Cryogenic optical resonator clocks are another type of 
clock used in tests of fundamental constants. Unlike 
atomic clocks, the rates of these clocks depend on the 
speed of light rather than Planck's constant. An invari- 
ant c means that the locally-measured (proper) rates of 
suitably configured resonator clocks will not vary with 
position, even if h varies. (Some resonator clocks are 
locked to an atomic frequency and would depend on h 
[l8|.) Therefore, a resonator clock measures proper time 
and can be used to determine /3(. 

Comparisons between atomic clocks (hydrogen masers) 
and cryogenic sapphire oscillators (CSO) have been per- 
formed at the Paris Observatory [J] and other laborato- 
ries (e.g., [iSf). Tobar et al. summarize that work and 
show that little measurable difference is observed as the 
earth moves in its elliptical orbit within the gravitational 
potential of the sun, within a limit of /3h maser — /3cso = 
—2.7(1.4) X 10^* for annual variations Q. The rela- 
tive rates of atomic and resonator clocks differ little with 
gravitational potential. 

Once a GPS satellite reaches orbit, the clock experi- 
ments [^ indicate that a resonator clock and an atomic 
clock on the satellite would advance at the same rela- 
tive rates. Therefore, the variations in clock corrections 
should apply equally to resonator clocks, and can be used 
to determine a limit on /?(. In other words, within a GPS 
orbit we can treat GPS atomic clocks as a surrogate for 
resonator clocks, so that Pt ~ /3/- 

In the clock comparison experiments, /3h masor — /3cso 
corresponds to l3f — (3t, and their limits were added to the 
limits on l3f to determine limits for /3t. Combining the 
maximum uncertainties, the limits on LPI violation for 
time dilation were found to be —0.0031 < j3t < 0.0037. 

To prove that h is invariant in the context of SME 
theory is probably not possible with existent experimen- 
tal evidence. However, in general relativity the locally- 
measured (proper) rest energy of a particle is indepen- 
dent of gravitational potential. Based on arguments sim- 
ilar to Nordtvedt's [20|, the proper energy emitted by a 
specific atomic transition should also be invariant. For 
example, the locally measured energy emitted by the 
transition between the two hyperfine levels of the ground 
state of Cs^'^'^ - the isotope used in atomic clocks - should 
be the same at all points and time in the universe. 

With a locally- invariant emission energy Eq, the 
proper frequency of emission fx of a given atomic tran- 
sition measured at any elevation x would then be given 
by 



Jx — ^o I I^X J 



(5) 



where dtx/dto is the ratio of the time rates at x and O, 
and fit is a dimensionless parameter that represents the 
extent of LPI violation for time dilation. If Pt = 0, the 
formula reduces to general relativity. 



where hx is the locally measured value oi h at x. If /i 
satisfies LPI, then hx would be constant and the proper 
frequency of an atomic clock would also be invariant. 
However, if hx varies with position then the proper fre- 
quency of an atomic clock would also vary. 



LPI violations for h can be written as 



hc,/ho = 1 + /?/, Af//c2 



(6) 



where h^ is the focahy measured value of h at reference 
point O, hj. is its locally measured value at x, and Ph is 
the parameter for Planck's constant. 

Let Jo be the proper frequency of an atomic clock when 
it is at O. Since Eq — hofo, it follows from Eq. (5) that 



Jx ^^ '^oJo/'^X' 



(7) 



Due to time dilation, the frequency fxo of the clock at 
X measured by an observer at O will be 



fxo = {dtx/dto){hofo/K) 



(8) 



A formula relating the three /3's to first order can be 
found by substituting Eqs. (1), (4) and (6) into (8), which 
yields /3h — Pt — Pf + 0(^^), where n = AU/c^. Second 
order relativistic effects in the earth's field are negligible, 
so higher order terms can be ignored, leaving 



3k~(3t-(3f. 



(9) 



The traditional interpretation of general relativity pre- 
dicts that I3t — (3f = 0, yielding /3/i = 0, which would 
signify that h is invariant. 

The estimated limits on /3t and /3f were used in Eq. (9) 
to determine LPI limits for Planck's constant. Combin- 
ing the limits on those two parameters and accounting for 
their signs, our final results, accurate to one significant 
digit, are summarized in Table L 



TABLE I: Estimated limits on LPI violation, based on this 
study unless otherwise noted. 



Type of LPI 



Limits found 



Basis 



Atomic clocks -0.0027 < /?/ < 0.0033 GPS data 
Clock comparisons /?/ - /3t = -2.7(1.4) x 10"* Ref. [4] 
Time dilation -0.0031 < fit < 0.0037 GPS & [4] 
Planck's constant |/3fej < 0.007 Eq. (9) 

To illustrate the effects of the filtering, we obtain 
\/3h\ < 0.01 if we use a 1-year filter but substitute an 
11-day moving average of /3/ (smaller than the 28-day 
lunar cycle) for the 51-day average and omit the 120-day 
filters. 

In conclusion, general relativity requires the local in- 
variance of the speed of light and emission energy. Within 



those constraints we find, based on experimental evidence 
from the GPS and clock comparison experiments, that 
h satisfies local position invariance to within a limit of 
0.007. The data also supports the gravitational time di- 
lation of general relativity for geometrodynamic clocks 
within a limit of 0.004. 

The work reported in this letter was carried out at 
California State University, Northridge. 
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